Quantum condensation in electron-hole bilayers with density imbalance 
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We study the two-dimensional spatially separated electron-hole system with density imbalance at 
absolute zero temperature. By means of the mean-field theory, we find that the Fulde-Ferrell state 
is fairly stabilized by the order parameter mixing effect. 
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In the imbalanced systems, where two spices of 
fermions with unequal densities interacts via interspices 
attraction, exotic quantum condensations, such as Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) and Sarma phase, 
are expected to appear. They are extensively investi- 
gated in various research areas such as heavy-fermion, 
organic superconductors in soHd state physics^, color su- 
perconductivity in high energy physics^ and ultracold 
fermionic atoms'^. The FFLO phase^fi- is characterized 
by the condensation of the pairs with a finite (nonzero) 
center-of-mass (COM) momentum. On the other hand, 
the Sarma (or breached pair) phase^'^ can be regarded 
as mixture of the quantum condensation of compos- 
ite bosons with zero COM momentum and degenerate 
fermion gas in the strong couphng regime. Such a rich- 
ness of quantum condensation phases in the density im- 
balanced systems present a great contrast to those in the 
usual density balanced ones, which show only a crossover 
between the Bardeen-Cooper-Shrieffer (BCS) condensa- 
tion in the week coupling regime and the Bose-Einstein 
condensation (BEC) in the strong coupling regime. 

Since 1960's, a great deal of efforts have been made to 
observe the quantum condensations of the electron-hole 
(e-h) pairs in strongly photo-excited semiconductors and 
semimetallic systems^. Among them, e-h bilayer systems 
are considered to be promising, since the radiative life- 
times of the statically separated electrons and holes are 
long enough to realize thermal equilibrium even at the 
extremely low temperatures. In fact, an indication of 
exciton BEC has already been found experimentally in 
coupled quantum wells (CQW) for the balanced case^. 
They also provide a stage for study on the exotic quan- 
tum condensations in the imbalanced systems, because 
e-h imbalance can be controlled by the gates and inter- 
layer bias voltagesi^. 

So far, most of the theoretical works on the quan- 
tum condensations of the e-h pairs have been focused 
on the balanced ones, and the studies on the imbalanced 
ones still lies far behind. Recently, Fieri et al. have in- 
vestigated the phase diagram of the e-h bilayer systems 
with density imbalance at absolute zero temperatureii. 
They solve the BCS self-consistent equations only for the 
Sarma phase, that is only for zero COM momentum of 
the pair, and identify its instability for the negative su- 
perfluid density as the appearance of the FFLO phase. 
Although they provide intriguing results, their inequiva- 
lent treatment on the FFLO and Sarma phases is insuf- 



ficient. 

In this paper, we investigate the phase diagrams of the 
e-h bilayer systems using the mean field approximation, 
which can take into account the Sarma and FFLO phases 
on an equal footing. As will be shown below, the resul- 
tant phase diagram is drastically changed: the Sarma 
phases are strongly suppressed and the FFLO phase oc- 
cupies relatively high density regions in the phase dia- 
gram. 

The Hamiltonian of our reduced model reads 



^fc.fc'4+q/2o'^-fc+q/2h'^-fc'+Q/2hCfe'+q/2o, (1) 



k,k' 

where c\,^ creates an electron (cr = e) or hole (cr = h) 
with wave vector k. The single particle energy of elec- 
trons/holes is expressed as ^fco- = A:^/2mcr — /io-, and 
the electron (hole) effective mass and chemical poten- 
tial are denoted by TOo (mh) and /io (/^h), respectively. 
We also introduce the interaction matrix element as 
Vk,k' = J v{r)e^^^~^ ^'^dr, where r is the in-plane 
coordinate, S is the area of the system, and v{r) is the 
interaction potential for the e-h attraction. In our Hamil- 
tonian, the electron-electron and hole-hole interactions 
are ignored, since we expect that their main role is the 
renormalization of the electron and hole bands. We also 
neglect the electrostatic (charging) energy accumulated 
among the layers and background charge, since such a 
Hartree contribution gives only a constant energy shift 
for the fixed electron and hole densities. 

Using the mean-field approximation, we obtain the 
thermodynamic potential at zero temperature as. 



= E ('^fe - + E ^^(fc) t^i^A,' A,(fc') 

k k,k' 

+ E [EiOi-Ei) + Ek^i-Ek)] , (2) 



'(77+)2 + A,(fe)2, = i?fe±ry„ 
-fe+q/2h)/2, Qix) is the Heaviside step func- 



where 

(Cfe+q/2o ± 

tion, and we introduce the order parameter Aq(k) with 
the COM and relative momenta q and k, respectively. 
For simplicity, we assume that the order parameters are 
real numbers. Because only the spatially homogeneous 
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order parameters are considered, the FFLO phase dis- 
cussed here is the "single-plane-wave" type, in which the 
wave function for the COM motion of the e-h pair is a 
plane wave. In the following, this type of FFLO phase 
is called simply FF phase, which is the original one pro- 
posed by Fulde and Ferrell^. 

The minimization conditions of the thermody- 
namic potential give the self-consistent equations, 
d^/dAq{k) — and d^/dq — 0, which determine the 
magnitude of the COM momentum q and the order 
parameter Aq(fe) for the given chemical potentials /Xg 
and /ih- The electron/hole numbers can also be calcu- 
lated via = —dil/dfia, and the extent of the im- 
balance is described by the polarization parameter a ~ 
{Ne — Nil) I {No + N-h). Solving these self-consistent equa- 
tions numerically and minimizing thermodynamic poten- 
tial Jl, we can identify (i) normal phase by Ag(fc) — 0, 
(ii) superfluid phase by Aq(fc) ^ 0, q = and a = 0, (iii) 
Sarma phase by Aq(fc) 7^ 0, g = and a 7^ 0, and (iv) 
FF phase by Aq{k) 7^ 0, g 7^ and a 7^ 0. 

We consider both the long- and short-range e-h inter- 
actions. As for the long-range interaction, we use the 
bare interlayer Coulomb potential, v{r) — /e^/r'^ + (P, 
where e is the electron charge, e is the background di- 
electric constant, and d is the interlayer distance. This 
choice of interaction is suitable for study of the usual 
CQW e-h systems, owing to the weak screening effects 
in the low-dimensional systems. In this case, we use the 
two-dimensional (2D) exciton Bohr radius oq = e/2e'^m^ 
and binding energy Eq — e^/eoo as the units of length 
and energy, respectively. We also set mh/nie — 4.3 and 
d = 2ao with the typical GaAs/AlGaAs CQW systems 
in mind. 

For the short-range interaction, on the other hand, we 
use the contact potential, v{r) = gS{r). To avoid the un- 
physical ultraviolet divergence, the coupling constant g 
should satisfy .g^"^ = S^^J2ki-^a + k^ /'2m,- — i5)^^ , where 
nil- = {m~^ + m^^)^^ is the reduced mass, Ea is the ex- 
citon binding energy, and (5 is a positive infinitesimal.^^ 
This contact potential is useful to consider the CQW e-h 
systems where the metallic gate structures are present 
nearby the quantum wells, since they screen the long- 
range Coulomb interaction. Such a consideration is also 
important because the Hamiltonian ([1]) becomes equiv- 
alent to that for the two-dimensional (2D) ultracold 
fermionic atom systemaA^-. Note that the derived self- 
consistent equations can be much simplified for the con- 
tact potential, because the order parameter has no k de- 
pendence. In this calculation, the units of length and 
energy are chosen as l/\/2m^Ea and Ea, respectively. 

Figure [1] shows the fj, — h phase diagram, where /i = 
{pe + fJ-h) /2 and h = — iJih) /2. For small /x (/i < 
—0.1), there is no electron and hole at /i = 0. This phase 
is referred to as the zero particle (ZP) phase. As \h\ 
increases at the fixed value of /i, the fully polarized (FP) 
normal phase with only electrons (a = -1-1), and that with 
only holes(a = —1) appear for positive and negative h, 
respectively. There is also another normal phase at large 
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FIG. 1: Phase diagram of CQW e-h system in the (/i, h) plane 
displaying superfluid(SF), FF, and three normal phases(PP, 
FP, ZP). Solid, dotted and dot-dashed lines denote the first- 
order transition, second-order transition and crossover, re- 
spectively. 



/i (/i > 0.6), which is the partially polarized (PP) normal 
phase characterized by \a\ 7^ 1. 

For > —0.1, the superfluid (SF) phase appears, 
which is the quantum condensation without polarization 
[a = 0) and has a region with a shape like a "needle" in 
the phase diagram. This "needle" has a positive gradient: 
the region of the SF phase is shifted to the larger positive 
side of h with increasing /x. This is because the chemical 
potential of electrons increase more rapidly than that of 
holes owing to the mass asymmetry mo < mh. 

One of exotic quantum condensation phases with a 7^ 
is the FF phase, which appears in the region around the 
SF phase. It is divided by the SF region [a = 0) into two, 
that are the electron-rich (a > 0) and hole-rich (a < 0) 
FF phases located at above and below the SF phase, re- 
spectively. The region of the hole-rich FF phase is wider 
than that of the electron-rich one, because the energy in- 
crease due to the reconstruction of the Fermi surfaces by 
the e-h pair condensation needs more energy for the elec- 
trons with light mass than for the holes with heavy mass. 
We can also see that the FF phase is much stabilized in 
the present system than in those with the short-range in- 
teraction, e.g. ultracold fermionic atom systems, which 
exhibit only a narrow FF phased. As discussed bellow, 
this stabilization can be explained by the order parame- 
ter mixing effect. 

Another exotic quantum condensation phase with a 7^ 
is the Sarma phase. This phase is expected to appear in 
the strong coupling region of the three-dimensional ultra- 
cold fermionic atom systems^ ! -'^^'-'^"^i-'^^'-'^^i-'^''''-'^^ . As shown in 
Fig. [1] however, it does not appear in the phase diagram 
for mh/me = 4.3, as well as the two-dimensional ultra- 
cold fermion systemai^iiSii^. Note that the Sarma phase 
is unstable for mh/rrie < 10 for both long- and short- 
range interactions. However, we find that this phase 
become thermodynamically stable only in the systems 
with the extremely large mass ratio, e.g. m/j/me = 30, 
and with the long-range interaction. This is because 
there is momentum dependence of the order parameter in 
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FIG. 2: Phase diagram of CQW e-h system in the (rs,a) 
plane. The abbreviations are the same with those used in 

Figm 

the case of the long-range interaction, as Forbes et al?^ 
claimed before. 

In order to compare our phase diagram with that 
obtained by Fieri et al^^ we redraw the phase di- 
agram on the (rs,a) plane, where the interaction 
strength is denoted by the dimensionless parameter — 
[tt (iVe -I- Nh) I2S\~^I'^. In Fig. H there are the SF phase 
on the line a = 0, the FF phase in high density regions, 
and "unstable" regions spread between them. In these 
"unstable" regions, there is no thermodynamically sta- 
ble and uniform (single-phase) solution. We can see that 
the most of the FF phase and all of the Sarma phase 
previously found are replaced by these regimes. This is 
presumably because only the positiveness of the super- 
fluid density is checked and the thermodynamic stability 
of the Sarma phase, i.e., the positive definiteness of the 
2x2 matrix dN„ / d^a' is not examined in Ref. |Tl|. 
If we ignore the effects of the background charge, it is 
expected that the phase separation takes place in the 
"unstable" region, as discussed in the ultracold fermionic 
atom systems^ii^ii^iiiii^ii^>^i2^. However, in the realistic 
CQW e-h system, it hardly occur since the inhomoge- 
neous charge distribution give rise to a huge electrostatic 
energy under the uniform background charge. 

The FF phase is spread in relatively high density re- 
gions and almost all of the regions corresponds to the 
normal phase in the previous study. The reason for this 
spreading is that Fieri et al. studied the self-consistent 
solutions only for the case of zero COM momentum of the 
pair. As is clarified below, this change is due to the mo- 
mentum dependence of the order parameter in the case 
of finite COM momentum of the pair. In Fig. [31 we show 
the fc-dependence of Aq(fc) for several choices of r^. and 
a. In the SF phase with a — 0, Ag(fc) shows isotropic s- 
wave behavior forming a circle on the {k^, ky) plane. As 
the polarization increases and the system enters the FF 
phase, the profile of Ag(fc) continuously deforms. The 
region with nonzero Aq(fc) in fc-space shrinks with in- 
creasing a and forms a shape of an arc. This clearly 
shows that the order parameter is not the simple s-wave 
but mixed wave in the FF phase, which is known as "par- 




FIG. 3: Momentum dependence of the order parameters. 
The upper panels are for n = 2.0 (rs — 0.5) and h = 1.24 
(a = 0.0), 1.19 (-0.04) and 0.79 (-0.30), from left to right. 
The lower panels are for /i = 0.5 {vs — 1.0) and h — 0.31 
{a = 0.0), 0.21 (-0.26) and 0.01(-0.61), from left to right. 
Solid and dashed lines are Fermi circle of electrons and holes, 
respectively. We choose the center-of-mass momentum as 
kx > and ky = 0. 
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FIG. 4: h dependence of the maximum values of Aq{k) for the 
long-range interaction A™'"' and A, for the short-range inter- 
action A™'''', which are nomalized by thir maximum values for 
each /i. (The subscript L and S indicate long- and short-range 
interactions, respectively.) The upper(lower) scales are for the 
short (long)-range interaction. 

ity mixing" or "order parameter mixing (0PM)" In 
the simple s-wave FF (sFF) state^i^^, the total energy 
of the system is decreased by the condensation but in- 
creased by the polarization, because the polarization be- 
tween sFF phase and normal phase are different for given 
chemical potentials. In contrast, in the mixed FF state, 
we have confirmed that there is no polarization difference 
between FF state and normal state, which indicates that 
there is no energy loss due to the polarization. 

As is shown in Fig. [U we also study the maximum 
values of Ag(fc) in the FF states for the long- and short- 
range interactions, which are referred to as A^^^ and 
^max^ respectively. They are normalized by their corre- 
sponding values in the SF phase. We clearly see that the 
behavior of Af^'^^ is quite different from those of A™^'^. 
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FIG. 5: /i-dependence of 5Q in the FF state. /iL(/is) in- 
dicates the mean chemical potential of the long(short)-range 
interaction. The upper and lower scales are for the short- and 
long-range interactions, respectively. 

In fact, A^'^^ in the SF regime is much larger than that 
in the FF state. As h decreases and the system enters 
the FF phase {h < 1.8 for fi = 6.0, h < 2.6 for ^ = 8.0), 
^max (decreases rapidly with increasing the polarization. 
In contrast, in the case of long-range interaction, A^'^'' 
is almost unchanged throughout the SF and FF phases 
due to the 0PM effects. 

Such behavior of the order parameters is directly re- 
lated to that of the thermodynamic potential. We define 
the thermodynamic potential Q measured from the cor- 
responding Q in the normal state as 5fl — (fi — r2o)/^^min, 
where fio and flmin are the thermodynamic potential in 
the normal state for given fj, and h and the minimum 
r2 — Oo for given /i, respectively. In Fig. [51 we show 
the /i-dependence of 6il for several choices of /i. For the 
short-range interaction, Sft almost vanishes in the FF 
phase, which indicates that the total energy gain in the 
sFF phase is very small due to the energy loss from the 
polarization difference and the small energy gain from the 
small Ag"^''. On the other hand, 5^1 for the long-range in- 
teraction shows significant values in the FF phase, which 
is about ten times as large as that for the short-range 
one. These results of the order parameters and the ther- 
modynamic potential clearly demonstrate the scenario 
discussed above. 



We also compare our results with those in the heavy 
fermion superconductors^^, where the simple SF and s-p 
mixed wave FF state are considered. At zero temper- 
ature, the ratio of the critical magnetic field in the s-p 
mixed FF state to that in the SF state has been obtained 
as 3.5. In our notations, this means that the ratio of the 
h range for fixed occupied by the FF phase to that oc- 
cupied by the SF phase is about 3.5. In our system with 
the long-range interaction, the ratio is 14 for fi = 2.0 in 
the electron- rich region (a > 0). In this sense, the FF 
state is much more stable in our system than in the heavy 
fermion superconductors. We can thus conclude that the 
CQW e-h system is a good candidate for the realization 
of FF state. 

Finally, we investigate the COM momentum q in the 
FF phase. We find that q/\kp — kp\ ~ 1.2 at the phase 
boundary between the SF and FF phase and q/\kp — 
kp\ ~ 1 at the boundary between the FF and normal 
phase. The similar behavior also appears in the case 
of the short-range interaction. Therefore, this behavior 
is independent of the range of the interaction and the 
0PM effects. On the other hand, the critical \a\ where 
the phase transition between the FF and normal phase 
occur is quite different between the cases of the long- and 
short-range interactions. The critical |a| is much larger 
for the long-range interaction than for the short-range 
one, which is due to the 0PM effects. 

In conclusion, we investigated the two-dimensional 
spatially separated electron-hole system with density im- 
balance at zero temperature. The results on the previ- 
ous work^^ is almost replaced by the present calculation, 
particularly in the high density region; The Sarma phase 
disappers and the FF phase are highly stabilized by the 
order parameter mixing effect in our system with the 
long-range interaction. We thus consider the CQW e-h 
system to be a good candidate for the realization of FF 
state. 
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